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Abstract 

Shuttle-assisted charge transfer is pivotal for the efficient energy transduction from the food-stuff 
electrons to protons in the respiratory chain of animal cells and bacteria. The respiratory chain 
consists of four metalloprotein Complexes (I-IV) embedded in the inner membrane of a mitochon- 
drion. Three of these complexes pump protons across the membrane, fuelled by the energy of 
food-stuff electrons. Despite extensive biochemical and biophysical studies, the physical mecha- 
nism of this proton pumping is still not well understood. Here we present a nanoelectromechanical 
model of the electron-driven proton pump related to the second loop of the respiratory chain, 
where a lipid-soluble ubiquinone molecule shuttles between the Complex I and Complex III, carry- 
ing two electrons and two protons. We show that the energy of electrons can be converted to the 
transmembrane proton potential gradient via the electrostatic interaction between electrons and 
protons on the shuttle. We find that the system can operate either as a proton pump, or, in the 
reverse regime, as an electron pump. For membranes with various viscosities, we demonstrate that 
the uphill proton current peaks near the body temperature T ~ 37° C. 
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Each mitochondrion has two functionally different membranes: the outer membrane, 
separating the contents of the mitochondrion from the cell cytoplasm, and the inner lipid 
membrane, embedded concentrically inside the outer envelope (see top right of Fig. 1). 
As Fig. 1 schematically shows, the inner membrane (in white), which is impermeable to 
charge and polar particles, divides the mitochondrial interior in two parts: (i) the central, 
negatively charged, compartment, referred to as the matrix (N, in yellow), and (ii) the 
intermembrane space (P, in blue), which is characterized by a positive charge. The electrical 
charge difference of the compartments results in an electrochemical proton gradient ~ —200 
mV jl, 2I, which is used to synthesize adenosine triphosphate (ATP) molecules, the main 
energy currency of the cell. This electrochemical potential is maintained by a set of proton 
oumps incorporated into the so-called Complexes I, III, and IV of the respiratory chain 
3, li, [sl, [7]. Proton pumps employ the energy of electrons, extracted from sugars and 
fatty acids after the digestion of food, for the uphill translocation of protons, from the 
negative (N) to the positive (P) side of the membrane. In the case of the respiratory chain 
in mitochondria, the N-side corresponds to the matrix, and the P-side is related to the 
intermembrane space (see Fig. 1). 

nn 

According to chemiosmotic theory |3j, |4j the pumping process in Complexes I and III 
is organized in three functional 'loops'. Each of these loops takes away some of the redox 
(electron) energy and converts it into energy of protons [8]. In the first loop, located in 
Complex I, two electrons are transferred from NADH (a reduced form of Nicotinamide 
Adenine Dinucleotide, NAD + ) to ubiquinone, pumping four protons in the process. 0, [h], 
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Our primary focus in this work is the second loop of the respiratory chain (see Fig. 1). 
There, a mobile molecule of ubiquinone, Q10 = C59H90O4, hereafter denoted by Q, takes two 
electrons from the NADH dehydrogenase enzyme (Complex I) and two protons from the 
N-side of the inner mitochondrial membrane, turning Q into a molecule of ubiquinol QH2. 
This molecular shuttle diffuses to the P-side of the inner membrane, where it donates both 
electrons to the cytochrome bci (Complex III). The oxidation of ubiquinol (and turning it 
back to ubiquinone) is accompanied by the pumping of two protons into the intermembrane 
space (P-side) against the proton electrochemical potential. Complex III transfers two elec- 
trons, acquired from QH 2 , to two molecules of cytochrome c, pumping in addition two more 
protons by means of the Q cycle mechanism (the third loop) 
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In this work we quantitatively analyze the proton pumping process in the second loop 
of the respiratory chain. Our goal is to show that the basic physical effect, related to the 
energetically uphill translocation of two protons per two electrons, transferred from NADH 
dehydrogenase, can be obtained within a simple nanoelectromechanical model, taking into 
account the electron-proton electrostatic coupling on the ubiquinone molecule (shuttle), 
diffusing inside of the inner membrane between Complexes I and III. This diffusion is gov- 
erned by the Langevin equation. In order to describe the process of loading/unloading the 
shuttle with electrons and protons, we employ a system of master equations, with position- 
dependent transition rates between the shuttle and electron/proton reservoirs. With these 
equations we analyze the time dependence of the pumping process, as well as the dependence 
of particle currents on both temperature and on electron/proton electrochemical gradients. 
There is a pool of ubiquinone/ubiquinol molecules in the inner mitochondrial membrane 
[5, Q], but we now only consider the contribution of a single ubiquinone molecule to the 
electron and proton translocation process. Because of this, the actual values of the electron 
and proton fluxes should be higher than the values calculated below. Note also that a similar 



proton pumping mechanism can work in photosynthetic systems 
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Model 



The mitochondrion can be considered as a nonequilibrium open system 2lj which has 
to be continuously supplied with electrons and protons. Complex I, characterized by a 
higher electron electrochemical potential, /ig = V e /2, serves as a Source of electrons (S), 
whereas Complex III, with an electrochemical potential /xd — — K/2, acts as a Drain (D), 
accepting electrons during the process of ubiquinol (QH 2 ) oxidation [8(. Note that we include 
the sign of the electron charge in the definition of the electrochemical potential (see also 
Methods). This means that a site with a higher electron energy is characterized by a more 
negative redox midpoint potential E m7 . The positive and negative sides of the membrane 
are characterized by the proton electrochemical potentials /xp = V p /2 and /xn = —V p /2. At 
physiological conditions, both electron (V e ) and proton (V p ) electrochemical gradients are 
positive: V e > 0, V p > 0. 

In our approach, the process of loading (from the source) and unloading (to the drain) 
the shuttle with electrons and with protons (from the N-side to the P-side of the inner 
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membrane) is governed by a set of master equations (see Methods) derived using methods of 



statistical mechanics and condensed matter theory |22|, 23], 24, 25]. The shuttle (ubiquinone 
molecule) has two electron-binding sites, l e and 2 e , with energies e\ e and e 2e , as well as 
two proton-binding sites, l p and 2 p , having energies €\ p and e 2p . To describe this situation 
quantum-mechanically we introduce 16 occupational states of the shuttle (see Methods). Our 
model for the nanomechanical proton pump can also work with just one electron and just 
one proton sites on the shuttle, in contrast to the kinetic models of pumps in cytochrome c 
oxidase 2ll. |25| . In its completely reduced form of ubiquinol QH 2 , the shuttle has a maximal 



load of two electrons and two protons, whereas in its oxidized quinone form (denoted by Q 
in Fig. 1) the shuttle is empty. 

The Coulomb interaction between charges located on the shuttle is given by the following 
parameters: the repulsion energy u e of two electrons, residing on sites l e and 2 e , and the 
repulsion energy u p of two protons located on sites l p and 2 p . The energy of the Coulomb 
attraction between an electron occupying the site a and a proton sitting on the site a' 
a, a' = 1, 2) is given by u ua i . In view of the spatial symmetry of the ubiquinone molecule 
2|, and for the sake of simplicity and without loss of generality, we describe all electrostatic 
interactions between electrons and protons on the shuttle by a single electrostatic energy 

= U , for any combinations of a and a'. 
We consider here the spatial projection of the shuttle motion on the straight line 
connecting the Q-binding sites of Complex I and Complex III. Along this line, the po- 
sition of the shuttle is characterized by a coordinate x. The Brownian motion [26J 
of the ubiquinone/ubiquinol molecule along this line is governed by an overdamped 
Langevin equation (see Eq. (3) in Methods) with drag coefficient ( and diffusion constant 
D = T/C (k B = 1). 

We find that for the optimal performance of the pump, all eigenenergies of the electron 
(E e = e ae ) and proton (E p = e ap ) sites on the shuttle (a = 1, 2) are approximately equal 
to one half of the charging energy Uo : E e = E p = Uq/2 > 0. This energy scale should 
be above the P-side proton electrochemical potential: Uq/2 > /ip = V p /2, in order to avoid 
discharging the proton battery through the shuttle. However, to allow loading the shuttle 
with the first electron, the energy, Uq/2, of the electron-binding sites on the shuttle should 
be near or below the potential, V e /2, of the source: Uq/2 < V e /2. This means that the 
charging energy Uq of the shuttle should fit the following transport window: V e > U > V p , 
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in order for the pump to work. 

The last cluster N2 of Complex I having a midpoint redox potential E m 7(N2) = 

—150 mV , can be considered as a source of electrons, and the Fe2S2-cluster of the complex III 
Q, Q], with a redox potential E m7 (Fe2S2) — 280 mV, is regarded as an electron drain. The 
difference between the redox potentials £ , m 7(Fe 2 S 2 ) and E m7 (N2) gives a rough estimate of 
the source-drain electron voltage buildup, V e ~ 430 meV. For the charging energy we select 
the value Uo = 450 meV, which slightly exceeds the electron potential difference V e ~ 400 
meV, and it is significantly higher than the proton voltage buildup V p ~ 250 meV [1,2]. We 
find that this choice, hampering the electron transport, results in almost equal electron and 
proton populations of the shuttle, and, thus, in a more efficient proton pumping across the 
inner membrane. 



B. Pumping mechanism 

The pumping mechanism can be illustrated by considering a single electron-proton pair 
(see Fig. 2). In this case the shuttle has one electron binding site (one electron potential well, 
shown in dark red in the upper part of each panel), and one proton site (one proton potential 
well, shown in grey in the lower part of each panel). The electron and proton charges of this 
pair almost completely compensate each other, and, because of this, they have practically no 
effect on the other electron-proton pair which can occupy the same ubiquinone molecule. In 
our example each electron or proton potential well has a single energy level (so-called "active" 
level, see black continuous lines crossing the wells). The dashed black line denotes the same 
electron or proton energy level, but shifted by the value of the electrostatic electron-proton 
energy Uq. This "virtual" level can be turned into the active level, and thus can be occupied, 
depending on the occupational state of the other well on the shuttle. To illustrate this, we 
introduce two possible energy scales, E e ±, E e2 , for an electron and two energies, E p i, E p2 , 
for a proton, residing on the shuttle. The higher electron (proton) energy E el (E pl ) points 
to the energy level of the electron (proton) for the case when the adjacent proton (electron) 
well is empty: E el = E e = U /2 (E pl = E p = U /2). The lower energy scale E e2 (E p2 ) 
gives the energy of an electron (proton) when the nearby proton (electron) potential well is 
occupied: E e2 = E e - U = -U /2 {E p2 = E p - U = -U /2). 

The shuttle, described by the position x, diffuses between the points ±xq. At x — +xq 
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(upper left corner of loop 2 in Fig. 1) the electron potential well is coupled to the source 
(S) and the proton well is coupled to the N-side of the membrane (Fig. 2a-2c). At x — —xq 
(bottom right of Fig. 1) the shuttle is close to the P-side of the membrane and to the electron 
drain (D) (Fig. 2d-2f). In Fig. 2 the sites populated with an electron or with a proton are 



shown as filled circles, w" 
transport theory 
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lereas the white circles show the empty sites. According to quantum 
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24i | . an electron (proton) can jump to the shuttle from the nearby 



electron (proton) reservoir if the electrochemical potential of the reservoir is comparable or 
higher than the active energy level of the electron (proton) on the shuttle. Unloading the 
shuttle with an electron (proton) to the nearby reservoir occurs in the opposite situation: 
when the electron (proton) energy level on the shuttle is higher than the potential of the 
electron (proton) reservoir. 

In the proton pumping regime, when V e > Uq > V p > 0, the electrochemical potential, 
Z^s — K/2, of the electron source is close to or higher than both possible electron energies E el 
and E e2 on the shuttle: /is > E e \ > E e2 . The electron potential of the drain, /id = —V e /2, 
is lower than the energy of the electron, regardless of the proton population of the shuttle: 
fii) < E e 2 < E e \. Therefore, the empty shuttle is always ready to accept an electron from the 
source, and the shuttle, populated with an electron, is ready to donate this electron to the 
drain. This means that the electron transport from the source to the drain is gated by the 
mechanical motion of the shuttle only, not by protons. In contrast, since /xn = — V p /2 < E p i, 
a proton can move to the shuttle from the N-side of the membrane only when the shuttle 
is populated with an electron, so that the proton energy is lower than the electrochemical 
potential of the N-side: E p2 < /in- The unloading of the proton cargo of the shuttle to the P- 
side of the membrane also occurs only when the electron previously escapes from the shuttle 
to the drain, so that the proton energy level increases above the electrochemical potential 
of the P-side: E p \ > /ip = V p /2. At this stage the electron energy is converted to the energy 
of the proton on the shuttle. Thus, besides the mechanical motion, the proton transfer from 
the matrix (N) to the intermembrane space (P) is essentially gated by the electrons. This 
gating prevents short-circuiting the proton transfer reaction, and, in addition, provides an 
efficient proton pumping across the inner mitochondrial membrane. 

The pumping cycle starts (see Fig. 2a) when the shuttle is near the N-side of the membrane 
and near the electron source S (x = +Xq). In Fig. 2b an electron jumps from the S-lead to 
the shuttle, since /is ~ E e i, thus lowering the proton energy on the shuttle from E p \ to the 
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level E p2 < /iN- In Fig. 2c a proton moves from the N-side to the shuttle, decreasing the 
electron energy to the value E e2 . Now the fully occupied neutral molecular shuttle diffuses 
to the opposite side of the membrane (from Fig. 2c to Fig. 2d). In Fig. 2d the shuttle is 
close to the P-side of the membrane and to the electron drain (D) ( X — — Xq ). In view of the 
relations E e2 > fir> and E p2 < /ip, the electron leaves the shuttle first, escaping to the drain 
reservoir (see Fig. 2e). Without the Coulomb attraction, the proton energy is increased 
above the level E p i, which is higher than the electrochemical potential fxp of the P-side. In 
Fig. 2f the proton jumps to the P-side, and after that the empty shuttle returns back to 
the N-side of the membrane and to the S-lead, thus completing the pumping cycle. As a 
result of this cycle, one electron is transferred downhill from the source (Fig. 2a) to the drain 
(Fig. 2f), and one proton is pumped uphill from the N (Fig. 2a) to the P-side (Fig. 2f) of 
the membrane. A similar mechanism describes the concerted transfer of two electrons and 
two protons (see Fig. 3 in the next section). 

C. Results 

To quantitatively describe the pumping process, we numerically solve the system of master 
equations together with the Langevin equation (see Eqs. (1,2) in Methods) for a parameter 
regime which provides a robust and efficient proton pump, and also corresponds to the 
ubiquinone molecular shuttle, Qi = C59H90O4. This molecule, which has the ability to bind 
two electrons and two protons, is widely distributed in the inner mitochondrial membrane 
of animal cells. 

In Fig. 3 we present the time evolution of the shuttling and pumping process for V e = 
400 meV, V p = 250 meV, and T = 298 K. The shuttle starts its motion at x = Xq (Fig. 3a) 
and after that diffuses between the membrane borders (shown by two dashed red lines at 
x = ±2 nm). On average, the shuttle crosses the membrane (the distance 2x = 4 nm) in a 
time At ~ 2.7 /is, which is very close to the characteristic diffusion time At = (2x ) 2 / (2D) 
for D = 3 • 10~ 12 m 2 /s. 

We find that, due to the symmetric configuration of the shuttle, both electron sites on 
the shuttle are populated and depopulated in concert: n e \(t) = n e 2(t). The same applies for 
the proton sites, n pl (t) = n p2 (t). Fig. 3b shows the time dependence of the total electron 
(n e = n e i + n e2 , continuous blue line) and proton (n p = n p \ + n p2 , dashed green line) 
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populations of the shuttle. At the beginning of the process (at t ~ and x ~ xq) the 
electron sites l e and 2 e are occupied simultaneously, with total population n e ~ 2. As 
a result of the concerted reduction of the shuttle with two electrons 18J, no intermediate 



semiquinone molecule, having a single unpaired electron, is formed. This conclusion is in 



agreement with the results of Ref. 



Ill ], where no semiquinone radical was detected in the 



reaction of the ubiquinone reduction by the tetranuclear cluster N2 of Complex I. With a 
tiny delay, the proton sites are also populated in concert. 

The almost neutral ubiquinol molecule QH 2 , loaded with approximately two electrons and 
two protons, diffuses and eventually reaches the P-side of the membrane (at t ~ 3 /xs, see 
Fig. 3a), where both electrons and both protons are transferred to the low-energy drain lead 
or to the high-energy P-side of the membrane, respectively (Fig. 3b). It follows from Fig. 3c, 
that changing the shuttle's populations in time is synchronized with the electron transfer 
from the source to the drain, described by the number N D (continuous blue line), as well as 
with the uphill proton translocation from the negative to the positive side of the membrane 
characterized by the proton number Np (dashed green line). The empty and neutral quinone 
molecule diffuses back, to the opposite side of the membrane (Fig. 3a), and the process starts 
again. Notice that, as a consequence of the stochastic nature of the Brownian motion, the 
shuttle is not populated and depopulated completely, and the proton population is a little 
bit smaller than the electron population of the shuttle (see Fig. 3b). Because of this, the 
shuttle acquires a tiny charge which can slightly increase the drag coefficient (. Here we 
ignore these changes. 

In Fig. 4 and 5 we show the average electron and proton particle currents, 
Id = (N d )/t r , I P = (N p )/t r , as functions of the electron (V e ) and proton (V p ) electro- 
chemical gradients at T = 298 K, diffusion coefficient D — 3- 10 -12 m 2 /s, and electrostatic 
energy Uq = 450 meV. The numbers of electrons, (iVo), transferred from the source (Com- 
plex I) to the drain (Complex III) in the time tr, and the number of protons, pumped 
from the N- to the P-side of the membrane in the same time t r , are numerically averaged 
here over several Brownian trajectories. Positive values of the electron current, Id > 0, 
correspond to the downhill (from source to drain) flow of electrons, whereas positive values 
of the proton current, Ip > 0, are related to the uphill translocation of the protons, from 
the N-side of the membrane (matrix) to the P-side (intermembrane space). At low values 
of V e and V p : V e < 200 meV, V p < 200 meV, both electron and proton currents are small 
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(Figs. 4 and 5). The pumping of protons against the potential gradient V p ~ 200 meV 
becomes possible at high enough electron potential gradients, V e > 300 meV (see the red 
area in Fig. 5), if V p < V e . At V e — 400 meV the system using a single ubiquinone shuttle 
can pump ~ 200 protons in one millisecond against the transmembrane potential V p ~ 250 
meV. This process is accompanied by a downhill electron flow (the red region in Fig. 4). 
Both electron and proton currents saturate at V e > 450 meV, if V p < 200 meV. It can be 
expected that the electron and proton currents across the membrane should grow almost 
linearly with increasing the number of the shuttles in the ubiquinone pool. 

For a sufficiently high electron driving force, V e > 450 meV, the electron current remains 
positive, Id > 0, even for a high proton voltage buildup, V p ~ 500 meV, whereas protons 
move downhill with I-p < 0. However, the proton component is able to reverse the electron 
flow at high proton gradients, V p > 300 meV, when the electron voltage is low enough, 
V e < 250 meV (blue area in Fig. 4). The reversibility of the electron flow in Complexes I 
and III in the presence of high proton transmembrane potentials has been observed in many 

nnnn 

respiratory and photosynthetic systems [10, lla, ll9|, l29(. In the reverse regime, the pump 



is able to translocate near 250 electrons per millisecond against the potential gradient of 
V e ~ 200 meV, if the proton potential gradient is V p ~ 400 meV. Figure 4 shows also that 
at sufficiently large electron driving forces (V e > 450 meV) the flow of electrons cannot be 
reversed (Id > 0) even by the very large transmembrane proton potential (V p ~ 500 meV). 

The temperature dependence of the electron, Id, and proton, ip, currents, at V e = 400 
meV and V p = 250 meV, is presented in Fig. 6 for V e = 400 meV, V p = 250 meV, and 
at three values of the drag coefficient C of the shuttle: Co/2 (blue continuous line, low 
viscosity), Co (green dashed line, intermediate viscosity), and 2Co (red dash-dotted line, high 
viscosity of the inner mitochondrial membrane). Here the drag coefficient Co = 1-37 nN-s/m 
corresponds (at T = 298 K) to the ubiquinone diffusion coefficient D = T/( = 3- 10~ 12 m 2 /s 



measured in Refs. 



27 



]. In order to gain further insight into the physical mechanism of 
the proton pump, we expand our plot to the two extreme regimes: very low and very high 
temperatures. It follows from Fig. 6a that for membranes with high viscosity (C = Co and 
C = 2Co) the electron transfer is accelerated and finally saturated with increasing temperature 
(Fig. 6a). In a warm environment, the shuttle performs more trips from the source to the 
drain and back transferring downhill more electrons. However, if the temperature is very 
high, and, especially, if the viscosity of the membrane is low (C = Co/2), electrons do not 
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have enough time to populate or depopulate the ubiquinone molecule. Thus, in spite of the 
more frequent trips, the shuttle carries less cargo, and the electron current decreases for very 
high temperatures and when ( = ( /2. This phenomenon is most conspicuous in the case of 
proton transfer (Fig. 6b). Loading (unloading) the shuttle with protons follows its loading 
(unloading) with electrons. Because of this, within the present model, the mobile shuttle 
moving in a warm medium has more chances to carry less protons, thus, slowing the proton 
transfer at high temperatures T > 400 K. As in the electron case, an initial increase of the 
proton current Ip with temperature reflects more frequent trips of the shuttle between the 
membrane sides. Note that, at any viscosity of the membrane, the uphill proton current has 
a smooth maximum near the physiological (or body) temperatures, T = 310 K = 37°C. 
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D. Conclusions 



We have shown that the electrostatic attraction between electrons and protons, travelling 
together on the molecular shuttle, can be responsible for the pumping of protons against 
the transmembrane electrochemical potential in the second loop of the respiratory chain 
(between Complexes I and III). The Brownian shuttle (a ubiquinone/ubiquinol molecule) 
has two electron-binding and two proton-binding sites. We have numerically solved a set 
of master equations, which quantitatively describes the process of loading and unloading 
the shuttle with electrons and protons, along with a stochastic Langevin equation for the 
shuttle position. We have found that, depending on the ratio between the electron potential 
gradient V e and the difference V p of proton electrochemical potentials across the membrane, 
the system can work either as a proton pump or, in the reverse regime, as an electron pump 
(Figs. 4 and 5). We have also studied the temperature dependence of the proton pumping 
for membranes with various viscosities and found that the uphill proton current is maximal 
near the body temperature T = 37°C (Fig. 6). 



E. Methods 



Occupational states of the shuttle. The occupational states of the shuttle is repre- 
sented by a set of basis functions, \e1e2 P1P2), where the occupations of the electron and 
proton sites e 1: e 2 ,pi,p 2 take the values or 1. The vacuum state, |1) = 1 00 00), corresponds 
to the empty shuttle (ubiquinone), and the shuttle loaded with two electrons and two protons 
(ubiquinole) is described by the state 1 16) = |11 11). These basis states, \p) (/i = 1, . . . , 16), 
and the spectrum of their eigenenergies, {E^}, are closely related to the electron-proton basis 



states and the energy spectrum introduced in Ref. 25] for the proton pump in cytochrome 
c oxidase. Here we do not consider the spin degrees of freedom of electrons and protons. 

Master equations describing the population of the shuttle. The probability of 
finding the system in the occupational state at the time t is determined by the diagonal 



element p^t) of the density matrix, averaged over the Fermi distributions of the e 
proton reservoirs. For these probabilities, we derive a system of master equations 



ectron and 
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25j,[30|: 
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with the relaxation matrix, 

D 2 
a=5 cr=l 

E Ea|W^(x)| 2 {|6^| 2 [l - + 

/3=7V <T=1 

l^/i^)}, (2) 

which depends on the position x of the shuttle and on the frequencies uj^ = E il — E v (h = 1). 
Here, T e aa = T e and T p 0a = T p are the energy-independent rates for electron and proton 
tunneling between the binding sites on the shuttle (a — 1,2) and the corresponding electron 
(a = S,D) and proton {(5 = N,P) reservoirs. Functions w a (x) and Wp(x) describe the 
dependence of the electron and proton transition rates on the position of the shuttle x. 
Also, (f and b° u represent matrix elements of the electron, a, and proton, b, annihilation 
operators, related to the electron or proton a— sites on the shuttle, and f a (E),Fp(E) are 
the Fermi distribution functions in the electron, /, and proton, F, reservoirs. The total 
probability is equal to one: X^=i p^ — 1- The populations of electron-binding sites, n e i, n e 2, 
and proton-binding sites, n p i,n p2 , of the shuttle can be expressed in terms of the functions 
p M . For example, for the population of the first electron site on the shuttle we obtain: 

n e i = p 2 + Pe + P7 + Pio + Pn + P12 + Pu + Pie- 
Using tools from statistical mechanics, we model the electron source and drain reservoirs 
by Fermi distributions over energy E with temperature T (ks = 1) and electrochemical po- 
tentials ps f° r the source and p^ for the drain: f a (E) = {exp[(E — p a )/T] + 1} _1 (a = S, D). 
The values of the chemical potentials, ps and pr>, provide the maximum energies of electrons 
in the last cluster of Complex I and in the first cluster of Complex III, respectively. Similarly, 
protons in the matrix (N-side) and in the intermembrane space (P-side) are described by the 
Fermi distributions, Fp(E) = {exp[(E — pp)/T\ + 1} _1 {(5 = N, P), with the corresponding 
proton electrochemical potentials p^ and pp. We choose symmetric configurations for the 
electrochemical potentials: ps = V e /2,p D = — V e /2, and p^ = —V p /2,p P = V p /2. The 
parameter V e is the difference in electrochemical potentials between the source (Complex 
I) and the drain (Complex III), and V p is the proton electrochemical gradient across the 
membrane. As mentioned before, the absolute value of the electron charge |e| is included 
to the definitions of electron and proton voltage drops V e and V p , which are measured in 
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meV. The positive values of the electron and proton voltages correspond to the case when 
/is > A*d and /i N < /-tp- 

Here we assume that the population and depopulation of the shuttle occurs via the tunnel- 
ing of electrons, characterized by the tunneling length A e . The Q— binding site of Complex 
I (the source of electrons) is located near the N-side of the membrane, with coordinate 
x = xq, and the Q— binding site of Complex III (the electron drain) is close to the P-side, 
with x = —xq. In this case we use the standard (in nanoelectromechanics, 23j) expressions 
describing the position-dependence of the source-shuttle and drain-shuttle electron-tunneling 
amplitudes: w$(x) = exp(— \x — xo\/\ e ), Wb(x) = exp(— |x + xo|/A e ). The region x > xq 
corresponds to the matrix (N-side), and the region with x < —Xq refers to the intermembrane 
space (P-side). Protons from both, N and P, aqueous sides of the inner membrane can move 
along proton-conducting pathways [17J and reach the Q-binding sites of Complexes I and III. 
Consequently, a proton exchange between the shuttle and the aqueous intermembrane space 
is possible when x > xq, and the proton transitions between the P-side and the shuttle take 
place when x < —xq- This situation is modelled here with the position-dependent proton 
amplitudes: W^(x) = [exp((xo — x)/\) + l] -1 , and Wp(x) = [exp((x + xq)/X p ) + l] -1 , 
characterized by the steepness X p . These forms quantify the above observations. 

The transition rate T e , between the electron-binding sites and the electron reservoirs, 
is assumed to be about the proton transition rate T p (between the aqueous sides of the 
membrane and the proton-binding sites on the shuttle): r e ~ T p ~ 5 ns _1 . For these rates 
the shuttle can be loaded (or unloaded) with electrons and protons in a time which is much 
shorter than the diffusion time scale, At ~ {2xo) 2 /(2D) ~ 2.7 fis, of the system. Here the 
distance 2xq = 4 nm is close to the membrane thickness. The electron tunneling length A e 
and the proton steepness parameter X p are chosen to be A e = X p = 0.25 nm. Protons can 
be delivered to the Q— docking sites by a system of proton-conducting pathways [l^ . We 
assume that the shuttle is loaded with N-side protons at x > (xo — X p ) , and unloaded in the 
P-side protons at x < (—x + X p ), where A p = 0.25 nm. The Q— docking sites are located 
on the opposite sides of the inner membrane and are separated by the distance 2xo- 

Stochastic equation describing the shuttle motion. Electron and proton transitions 
between the shuttle and the corresponding reservoirs have negligible effect on the mechanical 
motion. However, as follows from Eq. (2), electron and proton transition rates drastically 
depend on the spatial gap between the shuttle and the leads. The one-dimensional diffusion 
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of the Brownian shuttle along a line, connecting the quinone-binding site of Complex I and 
the quinone-binding site of Complex III, is described by the overdamped Langevin equation: 

dU(x) 

where ( is the drag coefficient of the shuttle in the lipid membrane. A zero-mean valued, 
(£) = 0, fluctuation force £ has Gaussian statistics with the correlation function: = 
2(T5(t — t'), which is proportional to the temperature T of the environment. A potential 
U(x) is responsible for the spatial confinement of the ubiquinone/ubiquinole molecule inside 
the inner mitochondrial membrane. The diffusion coefficient of ubiquinone, D = T/C, ~ 



28| at T = 298 K , corresponds to the drag coefficient 



3-10 m /s, measured in 
Co = 1-37 nN-s/m. We use the value of Co as a reference point (see Fig. 6). 
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FIG. 1: (Color online) Schematic diagram of the second loop of the respiratory chain. 

High-energy electrons are delivered to Complex I by NADH and forwarded to the ubiquinone 
molecule (Q) shuttling between Complexes I and III. Complex I serves as a source (S) of electrons, 
and Complex III plays the role of an electron drain (D). In addition to two electrons, the shuttle 
carries two protons taken from the matrix side (N) of the inner membrane at x = +xq. This 
fully loaded shuttle (the ubiquinol molecule QH2) diffuses freely to the P-side of the membrane, at 
x = —xq, and donates two electrons to Complex III. Thus, in this process, two protons translocate 
energetically uphill to the intermembrane space (P). 
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FIG. 2: (Color online) Energetic diagrams of the pumping process. The shuttle schemati- 
cally shown here carries a single electron and a single proton: S and D are the Source and Drain 
reservoirs for electrons (rust-orange), N and P are the Negative and Positive aqueous sides (blue) 
of the inner membrane. Of course, in reality the shuttle carries two electrons and two protons. 
However, it is simpler to explain the proton pumping process with only one of each. Potential 
wells for electrons (dark red, upper part of each panel) and for protons (grey, lower part) move 
with the shuttle from one side of the membrane (x = +xq) to the other side (x = —xq). In (a,b) 
the shuttle is near the N-side, loading an electron from the source (S). In (a) the energy level of 
the proton [see the lower part of (a)] is higher than the electrochemical potential of the N-side. 
In (b) the shuttle is populated with one electron. The effective energy of the proton is now lower 
than the potential of the N-side, and the shuttle is then loaded with the proton in (c). The fully 
loaded shuttle in (c) is still located near the N-side of the membrane (x = +xq). Both, electron 
and proton, lower their energies because of the electron-proton Coulomb attraction (vertical green 
arrows), (d) The diffusing shuttle carries its cargo to the P-side of the membrane at x = —xq. The 
energy level of the electron on the shuttle is now higher than the electron potential of the drain, 
and the electron translocates to the drain in (e). Without the electron-proton Coulomb attraction 
the energy level of the proton on the shuttle goes up, over the electrochemical potential of the 
P-side, and the shuttle gives up its proton in (f). Thus, the proton is pumped to the P-side of the 
membrane, in the lower part of (f), and the electron is in the drain. The shuttle is now empty and 
ready to diffuse back, to the N-side of the inner membrane, in order to start the pumping cycle 
again. 18 
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FIG. 3: (Color online) Time evolution of the pumping process, (a) Time dependence of 
the position x (in nm) (blue continuous curve) of the ubiquinone (ubiquinol) molecule Q (QH2), 
shuttling between the walls of the inner mitochondrial membrane located at x = ±xo (red dashed 
lines), where xq = 2 nm; (b) the total electron (n e = n e \ + n e 2, blue continuous line) and proton 
(rip = ripi + n P 2, green dashed line) populations of the shuttle versus time (in /xs); (c) the number 
of transferred electrons (Np>, blue continuous line) and the number of pumped protons (iVp, green 
dashed line) versus time at V e = 400 meV, V p = 250 meV, C = Coi an d at T = 298 K. Notice that 
the shuttle is loaded near the N-side of the membrane, at x ~ + xq, and unloaded at the P-side, at 
x ~ — xq. It follows from (c) that the process of shuttle unloading is accompanied by a stepwise 
increase of the number of electrons, Nd, transferred to the drain, and the number of protons, Np, 
pumped to the P-side of the membrane. 
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FIG. 4: (Color online) Electron current-voltage characteristics. The electron current Id 
(number of electrons transferred from the source to the drain per one millisecond) as a function 
of electron and proton electrochemical gradients, V e and V p , at Uq = 450 meV, £ = Co? an d at 
T = 298 K. This current is positive, Id > 0, at high enough electron voltage, V e > 300 meV, if the 
proton voltage is sufficiently small, V p < 250 meV. In the opposite case, when V p > 300 meV and 
V e < 250 meV, the pump works in the reverse regime, where the energy of the downhill-flowing 
protons (blue region in Fig. 5) is used for the uphill translocation of electrons against the potential 
V e (blue area in Fig. 4). 
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FIG. 5: (Color online) Proton current-voltage characteristics. The proton current Ip (num- 
ber of protons translocated from the negative, N, to the positive, P, side of the inner membrane 
per one millisecond) versus the electron, V e , and proton, V p , potential gradients, at Uq = 450 meV, 
C = Co> an d T = 298 K. The positive values of Ip (red area where V e > 300 meV, V p < 250 meV) 
correspond to the regime where the protons are pumped. Blue regions (V p > 300 meV, V e < 250 
meV) in Figs. 4 and 5 correspond to the electron pump regime, where protons move downhill 
(Ip < 0), thereby translocating electrons uphill, from the drain to the source. 



21 



500 




200 400 600 800 1000 



FIG. 6: (Color online) Effects of temperature and viscosity. Temperature dependence of (a) 
the electron current, Id, and (b) the proton current, Ip, in the proton pumping regime (V e = 400 
meV, V p = 250 meV) for three levels of membrane viscosity, determined by the drag coefficient 
(: (i) C = Co/2 (low viscosity); (ii) C = Co (intermediate viscosity); (hi) C = 2Co (relatively high 
viscosity). The drag coefficient Co = 1-37 nN-s/m is related (at T = 298 K) to the diffusion 



coefficient D = T/( = 3 • 10~^ m 2 /s 
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281 ] . In the region T < 250 K, both electron and proton 



currents increase almost linearly with temperature, because the shuttle performs more trips and 
translocates more particles as the temperature increases. At high temperatures, T > 600 K, the 
electron current (Fig. 5a) saturates and even slightly decreases (at low viscosity). The proton 
current (Fig. 5b) shows an even more pronounced maximum near T ~ 310 K (~ 37° C, body 
temperature), for all values of the membrane viscosity. Protons are loaded and unloaded on the 
shuttle after the electrons are loaded and unloaded. Because of this, at high T and low C) the 
mobile shuttle has not enough time to be loaded or unloaded with protons, carrying less protons 
across the membrane. 
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